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We derive computable upper bounds for the difference between an exact solution of the evolu- 
I tionary convection-diffusion problem and an approximation of this solution. The estimates are ob- 

tained by certain transformations of the integral identity that defines the generalized solution. These 
estimates depend on neither special properties of the exact solution nor its approximation, and in- 
• volve only global constants coming from embedding inequalities. The estimates are first derived for 

I functions in the corresponding energy space, and then possible extensions to classes of piecewise 

continuous approximations are discussed. 
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^ ■ 1 Introduction 
O 

I Let u{x, t) be the generalized solution of the initial-boundary value problem 



Ut = Cu iwQt-, ^(x, 0) = (1.1) 



where Qt '■= (0, T) x fi, is an open bounded domain in with the Lipschitz boundary dO., i2 is a 
linear uniformly elliptic operator, and 99 is a certain function that defines the initial condition. For this 
^ ■ class of evolutionary problems the existence and regularity theory has been deeply elaborated (see, e.g.. 

El 121). It is well known that if the boundary conditions on St '■= (0, T) x dQ are defined in a suitable 
way (e.g., u vanishes on St), then the generalized solution exists and is unique. Moreover, under some 
conditions imposed on the coefficients of C, 99, and dQ, one can establish the so-called first and second 
main inequalities that estimate different norms of u. 

In this paper, we focus on a different problem related to evolutionary models in mathematical physics. 
Assume that f is a given function (e.g., a numerical solution, or a solution of some simplified mathe- 
matical model associated with the same phenomenon) which we wish to compare with u. In general, 
the conditions that we impose on v are rather broad: we assume that it can be any function from the 
corresponding energy space. We do not suppose that v is subject to some extra regularity assumptions 
or that it must satisfy Galerkin orthogonality property with respect to some finite-dimensional space. 
Estimates of such a type are derived by purely functional methods without using specific properties of 
approximations and do not involve mesh dependent constants. For this reason, they are sometimes called 
functional type a posteriori estimates. 
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Finding computable and efficient estimates of a certain measure of u — v presents a mathematical 
problem of high importance for quantitative analysis of evolutionary problems in mathematical physics. 
We note that for linear and nonlinear elliptic problems (including stationary reaction-convection-diffusion 
problems and variational inequalities) this problem is well studied (cf. a systematic exposition in |7| and 
papers cited therein). In |61, it was shown that transformations of the integral identity that defines the 
generahzed solution of the parabolic heat equation lead to directly computable bounds for some weighted 
space-time norms of u — v. In [1], such estimates were obtained (and numerically tested) for a wider 
class of linear parabolic problems. Recently, in [5], similar estimates were derived for the evolutionary 
Stokes problem. 

In this paper, we consider the evolutionary convection-diffusion problem. We apply the same method 
as in [6 1, and show that for any v in the energy function space computable and guaranteed upper bounds 
of certain norms of u — v follow from the corresponding integral identity that defines the generalized 
solution. Moreover, we show that v may not exactly satisfy the initial boundary conditions (this matter 
is discussed in the last section of the paper) and it may be discontinuous in time. 

We consider the initial boundary value problem 

ut — divVn + a ■ Vu = f in Qt, (1.2) 
u{x,0) = (p{x) xen, (1.3) 
u = on St, (1.4) 

where Q is an open bounded domain in M^. We also assume that 



a = a (x) G L°° [n, Wj , diva G (Q) , diva < 0, (1.5) 
/GL2'i(Qr), ^eH^ifl). (1.6) 

We use standard definitions of function spaces associated with Qt', namely, L'^'^'{Qt) consists of func- 
tions from U{Qt), r >1, with the finite norm 

We denote by (Qt) the space {0,T;H^ (0)), and by (Qt) the subspace of (Qt) that 
consists of functions vanishing on St- The subspace V [Qt) consists of functions with the finite norm 



|t(;||y := ess sup ||w (•, t) ||^ + || Vi(;||2 
t6(o,r) 



Q 



T 



Functions in V (Qt) have bounded norms at cross-sections of the time-domain. The space 

V^'O (Qt) := {Qt) n C (O, T; (n)) 

is a subspace of V (Qt)- We denote by Vq^ (Qt) the subspace of V^'*^ that consists of functions vanish- 
ing on St ■ 

A function u G Vq'*^ {Qt) is called the generalized solution of the problem (ll.2l )- (ll.4l ) if it satisfies 
the following relation for all w G (Qt) 



Vu-Vwdxdt— / uwtdxdtA- I a-Vuwdxdt (1.7) 
Qt JQt JQt 

+ j {u {x,T) w {x,T) — u {x,0) w {x,0)) dx = / fwdxdt. 
Jn JQt 

The existence of the generalized solution follows from well known results (see, e.g., ||2]|3l). 
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2 Guaranteed Error Bounds for Space-conforming Approximations 

Letu € Hq {Qt)- In order to deduce a computable measure of u — ?;, we insert into the integral identity 
and represent it as follows: 

V {u — v) ■ Vw dx dt + / [a ■ V {u — v)) w dx dt — / {u — v)wtdxdt (2.1) 

Jqt Jqt 

^ V ' 

(*) 

+ / {{u{x,T) — V {x,T))w {x,T) — {u{x,G) — V {x,0))w {x,G)) dx 
Jn 

' V ' 

(**) 



: / {fw — Vv ■ \/w — vtw — {a ■ Vw) w) dx dt. 
Jqt 



We note that 



\w {x,t)[ 



t=T 



t=0 



dx 



\ll (^^ (^'*)) dxdt = j 



wwt dx dt. 



(2.2) 



Using (12.21) and setting w = u — v, from the two integrals (*) and (**) of (12.11 ) we find 



1 



n 



\u {x, t) — V (x, 



t=T f 

dx + \u {x,t) — V {x,t)\' 
*=o Jn 



t=T 



t=0 



dx 



4 / -v{x,t)f 
^ Jn 



t=T 



dx. 



t=o 



Hence, we arrive at the relation 



2 f 1 2 

|V(n — dxdt+ / [a ■ V {u — v j) {u — v) dx dt + - \\u — v\\q 
't JQt ^ 

(/ {u — v) — Vv ■ V {u — v) — vt {u — v) — {a ■ Vv) {u — v)) dx dt. 



t=T 

t=o 



Now, for w = w {x, t) we use the identity 

div [wa) = w diva + a ■ Vw, 



(2.3) 



(2.4) 



(2.5) 



which holds for t € [0, T] almost everywhere. We integrate it over Qt and take into account that w 
vanishes over its boundary St- We have 



(a • Vvu) w dxdt = j wa ■ Vw dx dt 
'Qt JQt 

Then, using (12.51 ) in the last term, we find 



div (wa) w dxdt. 



Qi 



Hence 



/ (a • Vw) w dxdt = — / (w^diva + wa ■ Vw) dx dt. 
Jqt Jqt 



{a ■Vw)w dxdt = —- j w'^ diva dxdt. 
't ^ JQt 



Using (12.61) . we rearrange (12.41 ) as follows 

f \V {u — v)\'^ dx dt / diva {u — v)'^ dx dt -\ — ||u — v| 

Jqt 2 2 



(2.6) 



(2.7) 



t=T 



f 1 2 

{f — vt — a -Vv) {u — v) dxdt — / Vv ■ V {u — v) dx dt + -\\u — v\\q^ 
Qt JQt 2 



i=0 



Let us define the error e = u — v. Since e = on 517, we have another relation for t G [0, T] almost 
everywhere 

/ {e diyy + y ■ V e) dx = / ey-jdx = 0, (2.8) 
Jn Jdn 

where 7 denote the unit outward normal vector to dQ, and the vector-valued function y = y{x, t) belongs 

to the space 

i?div(Qr) := [y G (0, T; R'^) ) , divy G (O, T; {^)) } . 

Then, using (12.81) in (12.71 ) we get 

/ (\Ve\^ -]-d\wa\e\^)dxdt + ]-\\e{-,T)\\l (2.9) 
Jqt ^ ^ ^ 

= / redxdt+ I {y - Vv) ■ Ve dx dt + - \\e {■,0)\\'^ , 
JQt JQt 2 

where 

r = r {v,y) := f — Vt — a ■ Vv + divy. (2.10) 
Henceforth, it is convenient to consider two cases; namely, diva = 0, and diva < 0. 

2.1 Case diva = 

If diva = (which is typical if convection is defined by a constant vector a) then the left hand side 
(LHS) of (12.91) contains the combined error norm 

lief :=l|Ve||J^ + ^||e(-,r)|l2 . (2.11) 



We use the Holder estimate 

'Qt 



re dx (it < llrllg^ llellg^ . (2.12) 



Since the following relation holds for t £ (0, T) almost everywhere 

||e||f^ < Cpn ||Ve||f^, 

we have 

llellg^ <C^J|Ve||Q^, (2.13) 
where Cfq is the constant in the Friedrichs inequahty. Therefore, 



'Qt 

Moreover, 



/ redxdt < Cp^ Mq^ IIVcHq^ < Cp^ MIq^ |||e|| 



/ {y - Vw) • Vedxdt < \\y - l|Ve||Q^ < \\y - ^v\\q^ |||e|| 

JQt 



(2.14) 



(2.15) 



Now (1221), (12141) . and (ITBl) imply 

|||e||p < ||y - Vv\\q^ M + WWq^ \\e\\ + ^ ||e (-,0)112 . 

Setting 

MQ^{v,y) := \\y-Vv\\Q^+CF^ \\A\q^. (2-16) 

and using some simple arithmetic-geometric inequalities, we deduce the simplest form of the guaranteed 
upper bound 

2 lllelll < Mqt {v, y) + (mI^ {v, y) + 2 ||e (•, ml ) ^'^ ■ (2-17) 
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It is easy to see that if the initial conditions are exactly satisfied, i.e. 

V (x, 0) = if (x) , 

then we get a simple form of the estimate 

M<MQ^{v,y), (2.18) 

which reflects the fact that the combined error norm is controlled by integral norms of discrepancies in 
the basic relations 

y — Vv = in Qt , 
divy + f — vt — a ■ \/v = in Qt ■ 

It is also possible to bound the term ||e (•, T) ||q alone. Using Young inequality in (12.141 ) we get 

/ redxdt< -7^^ Ml dt + / ^ ||Ve||^ dt, (2.19) 
Jqt Jo 2" -'0 2 

where a (t) is some positive bounded function, i.e. 

o<a{t)eL°° (o,r). 

Using Young inequality in (12.151) with some 

< /3(t) G (o,r), 

we get 

f {y-Vv)-Vedxdt< f l^\\y-Vv\\ldt+f ^WVeWldt. (2.20) 
Jqt jo -^P Jo ^ 

Then, dlD), (l2J9l) . and (IHOl) imply 
Choosing a = /3 = 1 we finally get 

\\e{-,T)f^<M\(v,y), {111) 

where 

M\ {v,y) := Cl^ ||r|||^ + ||y - Vv\\l^ + \\e (•, ml ■ 

2.2 Case diva < 

If -idiva = > 0, then we obtain estimates in terms of a different weighted norm 

l\4l-=\Ml^ + \Ml^ + \\e{;T)\\l. (2.23) 

By (IZ9l ) we find that 

\\VefQ^ + \\6e\\l^ + ^\\e{;T)\\l (2.24) 

<[ redxdt+ [ \\y-\7v\\^\\\7e\\^dt + l\\ei;0)\\l. 
JQt jo ^ 
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For reasons that will become clear later, we introduce a function X{x,t) with values in [0, 1] , and rewrite 
the first term on the right hand side (RHS) of (12.241 ) as follows 

/ redxdt= Xredxdt+ [1 — X)redxdt. 
Jqt JQt Jqt 



Then, for some < q (t) G L°° (0, T) we have 



Qt 



re dx dt < I ( - — 

4a 



< 



V4a 



\n il'^lln 



+ a\\6e\\l)dt+ / 11(1- A) r|| 

2 .T 

+ a\\6e\\l)dt + CF^ / \\{l-X)r\\^\\Ve\\^dt. 



(2.25) 



Using the Young inequality in the last term of (12.251) with some 

o<7(t) (o,r), 

we get 



/ redxdt < [ ( -— 
Jqt Jo ^4a 



A 



+ a ||5e||r, ] dt+ , , „ 
n ^ Jo ^ 27 



Hq_^)^||2 ^2||Ve||2 (2.26) 



Similarly, for some < /3 (t) G L°° (0, T), we have 



T 



T 



/3 2 /■ 1 2 

y-'^v\\^\\Ve\\^dt< I -\\Ve\\^dt+ tt^ l|y - Vv||f^ dt. 

JO ^ JO 



(2.27) 



By (I2.24H2.27I) . we conclude that 

(1 - f - |) II Ve||^ dt + (1 - a) Pe||^ dt + ^ ||e (•, T) 



(2.28) 



< 



1 

4a 



A 


2 ^ 




dt+ / 




n ^0 



'^^||(l-A)r||2 dt + 



j^\\y-Vv\\ldt + ^\\e{;0)\\l, 



where a(t), /3(t), and 7(t) are some arbitrary positive functions satisfying the conditions 

2-/3-7>0, l-a>0. (2.29) 
In particular, if we set a = /3 = 7 = 1 in (12.281 ). we can bound the term ||e (■, T) \\q alone as 



where 



A/?(A,i;,y) := 



|e(-,r)||^<Ar2(A,^,y), 



+ J|(l - A) r||^^ + ||y - VHI + lie (•, 0)11 



(2.30) 



If we set a = /3 = 7 = -, then we obtain the estimate in terms of the weighted space-time norm 

lllelll^ <2AA2(A,r;,y)-||e(-,0)||^. (2.31) 

Remark 1. The function A involved in these estimates can be useful if 6 attains small values at some 
points of Qt- In this case, ^ is a large penalty of the term r {v, y), and the estimate may be too pessimistic 
unless the value of r {v, y) is very small. In such a case, the function A (x, t) can be used to compensate 
this drawback. Indeed, we a priori know a and diva. Therefore, we can select A in such a way that j 



6 



remains of order 1. If 5 is sufficiently large every wfiere, tlien it may be useful to simply set A = 1 and 
use the simplified estimate 

£ (l - ^) II Ve||^ dt + «) \Ml dt + l lie (•, T) \\l (2.32) 



< 



1 

4a 



1 



2 r-T 

dt + 

n 



2/3 



^ ■■y-Vv\\ldt+^\\e{;0)\\l. 



Further, when the initial conditions are exactly satisfied, i.e. 

v{x,0) = (p{x), 

by choosing a = 1/2, and /? = 1, we have a simpler estimate 

1 



II III 2 ^ 

llellls < 



\y — Vu| 



Qt ■ 



(2.33) 



3 Nonconforming approximations 
3.1 Method 1 

Assume that we have obtained some nonconforming approximation v. The simplest approach to control 
its accuracy is to project it on a certain space of conforming approximations (e.g., by using averaging or 
other post-processing techniques). For elliptic type problems this method has been thoroughly discussed 
and tested in lH. Let P be some suitable mapping such that 

Then we can apply (12.281 ) (or (12.311 ) to this function. By the triangle inequality we obtain 

fu — v\lg < ll'U — -P^lll^ + lll^i — -P^^lll^ , (3.1) 

where the first term on the RHS is directly computable and represents the nonconformity error, whereas 
the second term can be estimated by the estimates derived in the previous section. It is worth noting that, 
in this approach, we do not exploit any specific structure of the underlying discretization method, and 
therefore, this approach is valid for any nonconforming approximation. 



3.2 Method 2 

We consider a special (but practically valuable) case of incremental approximations, in which the cor- 
responding function v is uniformly bounded and piecewise smooth (more precisely, it is smooth with 
respect to spatial and time variables on each time interval (t^, tk+i)) and may have jumps in time at 
points tfc. A; = 0, 1, ...A^, where to = and t^ = T. We show that the corresponding error estimate can 
be directly derived by a certain limit procedure applied to the main integral relation (I2.7I ). 
Indeed, from (12.71 ) we see that only the term 

— / vt{u — v) dx dt 
JQt 

involves the time derivative. To outline how the procedure acts in the context of time-nonconformity, we 
first consider the case where v has only one jump at a point t = r G (0, T). Assume that 

v{t — 0) = v{t + 0) = , ^ v'^ . 

On (0, T) we construct a sequence of functions v"^ (e is a small positive number) as follows: 

v\x,t) =v{x,t), t G (0,r - e) U (r,r), 

v^{x, t) = v{x, r - e) + - v') (*~^+"^ , t€iT-e,T) 



1 



It is easy to see that v'^{x, t) tends to v {x, t) in L (Qt) as e — > 0. Moreover, 

^ e( .X / Vt)(x,r), t e (0,r-e) 

' \ Vvix, r - e) + iVv+ - Vv') ^^±^, t G (r - e, r) 

and the space gradients also converge to Vv{x, r) in L^. 
We now apply the relation (12.71 ) to i;*^ : 



I V {u — v'') I (ix (it / diva (n — f ^) dxdt -\ ||n — v 



2.n^ ^ ' 2 



In 



(3.2) 



/ {f -vl-a-Vv^){u-v^) dxdt- I Vv^ -V {u - v^) dx dt + ]-\\u - v\\l^ 
Jqt JQt ^ 



Note that, as e 0, the LHS of (1X21 ) tends to 



f {V {u — v)f dx dt — - 1 diwa {u — i))^ dx dt + -\\u — v\\^ 
JQt ^ JQt 2 



t=T 



Moreover, 

dx dt, 



/ Vf^ • V (n - w') (ixdt ^ / Vv-V{u-v) 

JQt JQt 



IQt JQt 
and 

/ if - a ■ Vv") (u - v") dx dt / {f - a ■ Vv) {u - v) dx dt. 

JQt JQt 

It remains to consider the term involving the time derivative. We have 

/ / vludxdt—> j j vtudxdt+ / / v^udxdt (3.3) 
JO ^n Jo Jn ■ 



T—e 



r— e 



vtu dxdt + / ( H ) udxdt 

In Jo Jn Jr-e \ £ 

^ , , f v+ -V- r 

VfU dxdt + / / udxdt 

in Jo Jn f Jr-e 

vtudxdt+ I [v^ — u{t) dx, 
n Jo Jn 



since u (t) is continuous at t = r. Moreover, 



VtV^ dxdt= / VtV^ dxdt+ / v^v^ dx dt (3.4) 
'n Jo JnJo JnJr-e 

vtvdxdt+ I I dxdt 
n Jo Jn Jr-e 2 dt 

I I vtvdxdt + - I ({v^)^ — dx. 
JnJo 2 V / 



From (1331) and (134] ) it follows that 

(u — v*^) dxdt 



Qt 

vt {u — v) dxdt — / (f ^ — v") ru(r) 

Qt Jn ^ 2 / 

{)f (ti — ?)) — ^ (w^ — v^) [u (r) — v^) dx + - I (u^ — V )^ dx, 



Qt 



and we find that 



/ \Ve\^ dxdt-]- j d\wa{ef dxdt + ]-\\e{-,T)f^ (3.5) 
Jqt ^ JQt ^ 

■ I {f — vt — a ■ Vv) e dx dt + / Vv-Vedxdt 
Jqt JQt 

1 /" 1 

(-(;+ - V-) {u (r) -V-) dx + - / {v+ - v'f dx + - \\e {■, 0) ||^ 



m 2 7q 2 

< Cpn llrllg, II VellQ^ + ||y - Vv\\q^ || VeH^^ 



where f = r {v,y), and e = u — ?). Now recall that for n(r) — we can apply the majorant (12.221 ) or 
(12.30b depending on whether diva = or diva < 0, respectively. We proceed with both cases separately. 

3.2.1 Case diva = 

Note that u{t) — = e (■, r). Therefore, using (12.221) for u(r) — v~' we get 

h(^) -^~L = l|e(-,^)lln < {v,y) . (3.6) 
To simplify the notations, we introduce 

Cj = \\v+ - v-\\^ Mr {v, y) + \ \\v+ - v-\\l + \ \\e (•, ml ■ (3.7) 
Then, using (13.61) and (13.71 ) in (13.51) . and that diva = 0, we get 

llVellJ^ + \ \\e (•,r)||2 <[Cf,, Mq, + lly - ^v\\q, ) l|Ve||Q^ + Cj. (3.8) 
Since ||Ve||Q^ < |||e|||, we have 

llleli' < (C^o II^IIq, + lly - IIq^ ) lielll + Cj. (3.9) 
Using the definition (12.161 ) and some simple arithmetic-geometric inequalities, we finally obtain 

1 /2 

2\\e\\<MQ^{v,y)+[Ml^{v,y) + ^Cj) . (3.10) 

3.2.2 Case diva < 

In this case we use (12.301) for u{t) — v~, which gives 

||n(r)-?;-||^ = ||e(-,T)||f,<A/;(Ai,^,y), (3.11) 

where Ai {x,t) is a function with values in [0, 1] for t E (0, r). To simplify the notations, we now 
introduce 

'^J = \y - v'L^r {Xi,v,y) + \ - v-\\^ + \ ||e(-,0)||^ . (3.12) 
Proceeding in the same way as in the conforming case, we finally obtain 



£{1-^-1) \\Ve\\l dt + «) \ml dt + \\\e (., r)||^ (3.13) 



<Vj + 



' dt+ -^\\il-X)r\\ldt+ \\y-Vv\\ldt. 
Q Jo 27 Jq 2(j 



4a 



A. 



If we set a = /3 = 7 = ^, then we get a simpler estimate 

III - III 2 



e\li < 2Vj + 2Mi (A, v,y)-2 \\e (•, 0) ||^ . (3.14) 



Remark 2. If, instead of (I2.23I ). we introduce a slightly different weighted norm 



elllf := \Ne\\l^ + \\6e\\l^ + ^\\e{;T)\ 



(3.15) 



then, using (13.111) and (13.121) in (1331) . and that || Ve 



Qt 



< III e III ^, we get 



ellll < {Cf, 



r 



Qt 



+ \\y — \7v 



Qt 



(3.16) 



which easily leads to the following estimate similar to (13.101 ) 



2 lllelll^ < Mq^ {v, y) + [mI^ (z), y) + 4Pj) 



1/2 



(3.17) 



We see that the estimates in both the cases (through (13.71 ) and (13.121 )) involve penalty terms depending 
on the jump \\v^ — v~ ||q. If the latter quantity is small, then the overall value of the majorant does not 
essentially increase. If we wish to introduce more time-discontinuity points then we can easily extend 
these estimates using the techniques discussed to compute ||e (•, T) ||q in Subsections 12. 1 1 and 12. 2 1 

Remark 3. It may be very convenient to use approximations discontinuous in time if the spatial dis- 
cretizations are reconstructed during the process of time integration. In this case, at certain time moment 
r, we need to change the structure of the finite-dimensional space which is used for approximation of 
the spatial components of the solution. Then, we may have difficulties in conforming continuation of 
the approximate solution. With the help of nonconforming extensions (which are technically simple) 
we can easily obtain f (r + 0) by interpolating f (r — 0). Jumps of discontinuities that arises in such a 
procedure can be taken into account due to the penalty terms in the above estimates. 
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